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Abstract. We discuss possible sources of biases in the determination of the 
initial mass function (IMF) introduced by the binning of the data, the uncer- 
' tainty in the determinations of masses, and the existence of unresolved multiple 

systems. Those three effects tend to produce IMFs that are flatter than the real 
one. We analyze the importance of each effect and suggest techniques that mini- 
mize or eliminate the biases. We also report the detection of the first astrometric 
binary system composed of two very-early O-type stars, HD 93129 A. 
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S^. ■ 1. Introduction 
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In order to obtain the mass function of a stellar population from photomet- 
ric data one starts by placing the stars in a color-magnitude or in a theoretical 
(temperature-luminosity) HR diagram along with the evolutionary tracks and 
the corresponding isochrones. For a simple population one can then find the 
appropriate isochrone and obtain the masses for each star. For a complex star 
formation history, one has to obtain the full transformation from temperature 
and luminosity (or color and magnitude) to mass and age. The number of stars 
as a function of mass (corrected for age effects if necessary) can then be used to 
obtain the initial mass function (or IMF). 

A number of problems and biases can arise along the way to obtaining the 
true IMF due to the oversimplification of the assumptions (e.g. using isochrone 
fitting for a population with an age spread), the existing intrinsic degeneracies 
(e.g. the complicated topology of evolutionary tracks, metallicity effects, and 
rotation can make two stars of different masses and ages have the same temper- 
ature and luminosity), as well as for other reasons. In this work we will analyze 
three sources of systematic effects: the numerical bias introduced by the use of 
constant-size bins for the fitting of the IMF, the "mass diffusion" from low to 
high masses due to photometric uncertainties, and the existence of unresolved 
multiples. Those three effects go in the same direction of making the measured 
IMF flatter than the real one and can affect different samples to different de- 
grees, thus introducing the possibility of yielding a dispersion of measured IMF 
values where only a single one exists in reality. 
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2. Binning biases 

The first type of bias we will discuss is purely numerical and affects not 
only IMF determinations but the fitting of any function to binned data that 
follows Poisson or multinomial statistics ( Wheaton et al 1 ll995l : iLucvl fcoom and 



was analyzed bv IMaiz Apellaniz <fe Ube"dal (2005). Suppose we have measured 



the masses (m) for a set of stars to which we want to fit a power law of the form: 

^- = A-m\ (1) 
dm 

where dn is the number of stars with mass in the interval m to m + dm. In- 
tegrating both sides of the equation and taking logarithms we arrive at the 
expression: 



A 



10 



7 + 1 



Ami \ 7+1 ( Am, : \ 7+1 



(2) 



where re, is the mass at the center of an interval of width Am, that contains Ni 
stars. Equation |21 is the expression that can be used to derive the IMF slope, 7, 
by measuring the number of stars in each bin and fitting the data, which can 
be done by minimizing a x 2 statistic and deriving the associated uncertainties. 
Since Ni follows a binomial distribution (with N being the total number of stars 
summed over all bins), the associated weight for bin i for the x 2 fit is given by: 

NiN 

Wl ~ (JV-JViXlogioe) 2 ' (3) 

Two warnings should be given here. The first one is that for two bins i and 
j, Ni and Nj are cor related beca use they are both part of a joint multinomial 
distribution (see e.g. lLucvll200 0). The second and most important one is that 
the value of Ni that should be strictly used i n Eq. El is the one cal culated from 



the fit, not the one measured from the data ( Wheaton et al.lll995h . Otherwise, 
one runs into the possibility of introducing biases in the measurement of 7 unless 
some precautions are taken because of the differences between the real and the 
assumed weights assigned to bins with a small number of stars in them. 

Using the fit Ni instead of the data Ni for the weights requires iterating and, 
therefore, complicates the IMF calculation. An alternative strategy was analyzed 
bv lMaiz Apellaniz Ubedal ((20051 ): rather than trying to find the right weights 
for each bin, one can select the bins in such a way that they all have similar 
weights and, therefore, the obtained value of 7 is nearly independent of the 
weights themselves. In this manner, the weights derived from the data should 
give very similar results to the weights derived from the fit, hence eliminating 
the need for an iterative procedure. This strategy can be implemented by using 
bins with object equipartition, i.e. selecting the bins in such a way that they all 
have the same number of stars in them. A graphical example of the difference 
between constant bin-size (the standard approach to fitting binned data) and 
variable bin-size with object equipartition can be seen in Fig. ^ 

In order to test whether variable bins provide a significant improvement over 
constant bins, IMaiz Apellaniz Sz Ubedal (|2005h designed a series of Montecarlo 
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Figure 1. Comparison between the data and the fitted functions for one of 
the realizations with 30 stars and 10 bins for the two numerical experiments 
described in the text. The left panel corresponds to the first experiment 
(constant bin size, fixed lower and upper mass limits) while the right panel 
corresponds to the second experiment (variable bin size with object equipar- 
tition, lower and upper mass limits determined from the data). Note the 
differences in the size of the error bars between the two experiments and that 
the left panel plot includes four bins with zero stars. 



Table 1. Normalized biases for the first experiment (constant bin size) for 
3, 5, 10, 30, and 50 bins. 
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Table 2. Normalized biases for the second experiment (bin equipartition) 
for 3, 5, 10, 30, and 50 bins. 
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Figure 2. (left) Histogram with the distribution of (7^ + 2.35) juh for the 
1000 realizations of the first experiment with 300 stars and 50 bins. A Gaus- 
sian distribution with mean b — 2.200 and dispersion of 1.0 is also plotted for 
comparison. The vertical lines mark the position of and of b. (right) Bias as 
a function of -/Vi iin i n for the first experiment. Note that Ni tm [ n can be smaller 
than 1 because it is a property derived from the parent distribution. 
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Figure 3. (left) Histogram with the distribution of (7^ + 2.35)/cfc for the 
1000 realizations of the second experiment with 300 stars and 50 bins. A 
Gaussian distribution with mean b — 0.090 and dispersion of 1.0 is also plotted 
for comparison. The vertical lines mark the position of and of b. (right) 
Bias as a function of Ni for the second experiment. Note that the vertical 
scale for the plot is 1/10 of that of the right panel in Fig. [3 
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numerical experiments in which they also analyzed the importance of the choice 
of the lower and upper mass limits. The reader is referred to that article for 
details; here we provide a summary in Tables ^ an d HI and in Figs. |2 and EH of 
two of the experiments. In the first experiment, they generated 1000 realizations 
from a distribution with Salpeter (7 = —2.35) slope using either N = 30, 100, 
300, or 1000 stars. Each realization was then binned into 3, 5, 10, 30, and 50 
constant-size bins and a power law was fitted using \ 2 minimization. Finally, for 
the 1000 realizations of each N + number of bins combinations the normalized 
bias b was computed using the definition: 

, 1000 . ok 

1000 fc ^ u k K ' 

where is the uncertainty in 7^ derived from \ 2 minimization and the k index 
is used to denote the realization number, b is a sensible choice to judge the 
existence of biases. If |6| <C 1, then the fitting method will be unbiased because 
it will yield values that will be larger than the real one on ~50% of the occasions 
and smaller on another ~ 50%. If, on the other hand, \b\ ~ 1 or larger, a 
significant bias will exist. 

As seen in Table ^ significant biases exist for most N + number of bins 
combinations when using constant bins. In the right panel of Fig. [21 we see that 
b is a strong function of iVj im i n , the mean iVj in the bin with the lowest number 
of counts (which for constant-size bins and a Salpeter power law will be the 
rightmost one). Note that at least 10 stars in the rightmost bin are required 
for this method to yield small biases, thus making its use impractical for most 
applications. Also note that the bias always points in the direction of making 
the measured IMF flatter than what it really is and that, for a fixed number of 
bins, the effect is larger when there are fewer stars. 

The second experiment (third in the numbering of lMaiz Apellaniz &: Ubedal 

2005) was a repetition of the first one using bin equipartition instead of constant 
bins. As a comparison of Tables ^ and [2 shows, the biases for the second exper- 
iment are much smaller (even for the extreme case of dividing 30 stars into 30 
one-star bins, b is only 0.264). The effect is also perceptible in the right panel 
of Fig. |SJ where the vertical scale is a factor of 10 smaller than in the equivalent 
panel of Fig. |2J Furthermore, the left panel of Fig [2 shows that the distribution 
of (7^ + 2.35)/o"fc is well approximated by a Gaussian with a dispersion of 1.0, 
indicating that the random uncertainties derived from \ 2 minimization can also 

be trusted. 

W e conclude that the bin equipartition strategy proposed bv lMafz Apellaniz Ubedal 
( 2005) for the fitting of power laws with Salpeter slopes yields results that (a) 



are nearly bias-free and (b) produce correct uncertainty estimates. On the other 
hand, the standard uniform-size binning introduces biases that are dependent on 
the number of stars per bin. The power of the equipartition technique extends 
to small samples, since it is possible to obtain accurate values with reasonable 
precisions for the IMF slope even when as few as 30 stars are available for anal- 
ysis. 

We would also like to point out that, given the purely numerical nature 
of the analysis, these results could be extended to other similar problems. For 
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example, the mass function for young stellar clusters ca n be rather well ap - 
proximated by a power law with a slope of —2.0 (see e.g. iFall h Zhandl200lT ). 
which is quite close to —2.35, so the same type of biases should be present there 
as well. In general, we recommend that biases be evaluated for any function 
fitted to binned data through % 2 minimization by means of specific numerical 
experiments analogous to the ones in this article. 



3. Mass diffusion 

The second type of bias we will discuss is caused by photometric uncer- 
tainties and detection limits. When one tries to use photometric data to derive 
statistical properties of a stellar population, one finds that corrections for the 
undetected stars must be included in the calculation because it is easier to de- 
tect bright stars than dim ones. Such an incompleteness correc tion is usually 
handled by c rowded-field ph otometry packages such as DAOphot (|Stetsonlll987l ) 
or HSTphot ( Do lphinll2000h by doing experiments in which the code is run with 



artificial stars added and the percentage of recovered objects as a function of 
magnitude and color is calculated. The correction is then applied to the observed 
stellar statistics. 

Obviously, ignoring an incompleteness correction can result in a large bias in 
the measurement of the IMF and this well-known fact is taken into consideration 
in modern works on the subject. However, a related bias which is more subtle 
is not always taken into account. Suppose that we observe several times a 
star that has a real magnitude m. Due to Poisson, detector, and background 
noise, in some of our observations we will measure a magnitude m' > m and 
in others we will measure m! < m. If our detector is well calibrated, the first 
circumstance will happen 50% of the time and the second one the remaining 
50%. From our analysis of the detector we should be capable of estimating from 
a single measurement of the star an uncertainty o m in such a way that m! ± a m 
behaves in an approximately Gaussian way, e.g. the single-measurement values 
m' will be within m — a m and m + a m for approximately 2/3 of the sample 
and outside (m — 2a m ,m + 2a m ) for approximately 5% of the sample. Now, 
for most of the stellar mass range the IMF has a negative slope, meaning that 
there are more low-mass (dim) stars than high-mass (bright) ones. Therefore, if 
we measure a star to have magnitude m', there should be a higher probability 
that its real magnitude m is dimmer than m' than that it is brighter (i.e. there 
are more dim stars disguised as bright ones at a given measured magnitude 
than bright stars disguised as dim ones) because the underlying real luminosity 
distribution provides more dim stars to start with. We will call this effect mass 
diffusion because it acts in a manner analogous to a diffusion process, smoothing 
a gradient by shifting objects from where they are more abundant to where they 
are more scarce. 

Computing the correction required to eliminate mass-diffusion effects from 
an IMF is not straightforward because two intermediate steps are required. First, 
one has to translate uncertainties in the measured magnitudes and colors into 
uncertainties in temperature and luminosity (or bolometric magnitudes). Sec- 
ond, the uncertainties in the theoretical HR diagram have to be converted into 
uncertainties in mass. The first step involves applying extinction and bolometric 
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Figure 4. Converting from photometry-derived temperatures and luminosi- 
ties to masses. This plot shows the temperatures and luminosities for two stars 
in NGC 4214 derived from HST/WFPC2 F170W+F336W+F555W+F814W 
photometry and th eir associated uncertainty ellipses. Also shown are evolu- 
tionary tracks from iMaeder fc Mevneti l|2001f) . 
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corrections and taking into consideration the correlations between them and the 
measured magnitudes. An example of such effects is shown in Fig. where we 
have plotted the temperatures and bolometric magnitudes of two stars in NGC 
4214 derived from multiband HST/WFPC2 stellar photometry (|Ubeda et alJ 
2005). A strong correlation is observed between temperatures and bolometric 
magnitudes: most of this correlation is caused by the strong temperature depen- 
dence of the bolom etric correction. The data plotted in Fig. 0] was calculated 
using CHORIZOS ( Maiz Ad ellanid 1200 41. a code specifically designed for the 
task of transforming from measured magnitudes to physical properties such as 
temperature, age, or extinction. Note that the lower uncertainty ellipse is larger 
than the upper one: most of this effect is caused by the higher uncertainties in 
the measured magnitudes of dim stars compared to bright ones. 

The second step (calculating the uncertainties for the masses) can be achieved 
by producing a point-to-point coordinate conversion between temperature + lu- 
minosity and mass + age (with the caveats about one-to-one correspondence 
previously mentioned), generating a distribution of temperatures and luminosi- 
ties for each star according to its uncertainty ellipse, transforming those values 
into masses and ages using the conversion above, and calculating the mean and 
standard deviation of the derived mass distribution for each star. Those values 
can then be used as the mass (-Mi) and its uncertainty (om^) for each star. 

Once the masses and corresponding uncertainties have been computed for 
each star, one can calculate the correction due to mass diffusion in the following 
way. First, a function gm[M) is computed from the values for the individual 
stars (Mj, 0mJ by fitting a simple function such as a parabola. Then, one 
generates a series of realizations of the IMF with different values of the slope 
(which we will call 7 rea i) which are then smoothed using a Gaussian kernel with 
ctm{M). The output is fitted using x 2 minimization and the corresponding fitted 
value of the slope (which we will call Tat) is obtained. Finally, as shown in Fig. [51 
a polynomial is fit to 7fi t — 7 re ai as a function of 7fit which is the correction that 
needs to be applied. Note that since <ta/(M) is data-dependent, an individual 
correction has to be applied to each specific observation. 

We show in Fig. |S] the correction for the specific case of the NGC 4214 
data previously mentioned. The magnitude of the correction can be taken to be 
typical for HST photometry beyond the Magellanic Clouds obtained with a few 
orbits of exposure time. The bias produced by not applying the correction is in 
the same direction as the one caused by constant-size bins: the measured IMF 
appears to be flatter than the real one. 

We recommend that the correction described here be applied to the calcu- 
lation of IMFs in general. However, we should point out that, ideally, one would 
like the correction to be as small as possible. One (obvious) way to achieve 
this is to obtain photometry with better S/N ratio. Another one is to use as 
many filters as possible to adequately characterize the temperature (and possi- 
bly gravity and metallicity) of each star and to adequately correct for extinction 
and then to process the data using a code like CHORIZOS, instead on rely- 
ing on conversions from single-color + magnitude diagrams to temperature + 
luminosity equivalents. 
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Figure 5. Correcting mass-diffusion biases. The plot shows the difference 
between the fitted and the real slopes of the IMF as a function of the fitted 
value for the artificial IMF realizations described in the text. The individual 
data points are the results of the fit to the realizations and the continuous 
line is a polynomial fit. The <7m(M) and the total number of stars used are 

those corresponding to the NGC 4214 data of lUbeda et all l)2005(K 
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4. Multiplicity 

The third type of bias that will be discussed here is that caused by mul- 
tiplicity. It has been known for a long t ime that a large fraction o f stars are 
located in multiple systems. For example, iKouwenhoven et al.l (|2005h measured 
that at least 61% of the stars in the Scorpius-Centaurus OB association are in 
a multiple system (the number could actually be higher due to incompleteness 
and selection effects). Interestingly, the single-star fraction decreases for early 
spectral types: those authors measured that 41% of the systems in which the 
primary is a B4-B9 are single but the fraction decreases to 12% if the primary is a 
B1-B3. The large multiplicity at the hi gh-mass end of the stellar mass spectrum 
is not a new result: iMason et al.l ( 1998h measured that at least 75% of the stellar 



systems with O stars in clusters or associations are multiple. Furthermore, those 
authors recognized that with the current instrumentation capabilities there is 
still enough discovery space between visual and spectroscopic binaries to al- 
low for basically all O stars in clusters and associations to be part of multiple 
systems. 

The existence of unresolved binaries artificially flattens the IMF due to a 
combination of two effects. First, an unresolved binary of any mass is shifted 
from a lower mass to a higher mass in a mass histogram. Second, the flattening 
should be further enhanced by the apparent increase of multiplicity with mass 
in the range between a few and several tens of solar masses , which is the range 
for which Salpeter slopes are reported by most authors ((Chabrierl 120031 and 
references therein). 

The flattening due to unresolved binaries has been estimated by iKroupal 

1 200lh to produce a change in 7 between 0.0 and 1.3 for low-mass stars and 



brown dwarfs. Given the large multiplicity fractions observed for OB stars, the 
effect must also be significant for intermediate- and high-mass stars. At the 
highest end of the stellar mass spectrum, the problem of the calculation of the 
IMF slope is coupled with another one: is there an upper limit for stellar masses 
or is the highest mass in a cluster simply determined by statistical sampling in a 
quasi-Salpeter power law that extends to infinite masses? A few years ago there 
was a large discrepancy between the highest mass measured from orbital motion 
(the most reliable method of measuring masses), which yielded values around 
60 Mq, and the masses measured from photometry and spectral classification in 
R136 bv lMassev &: Hunterl (1998). Those authors measured stellar masses in the 



range 120-150 M Q and claimed that the data were compatible with a Salpeter 
IMF that extended beyond there. That gap has been recently narrowed in both 
directions. On the one hand, the 60 M barrier for spectroscopic-binary masses 
has been broken and the current heavyweig ht champion, WR 20a, lies around 80 
M (|Rauw et al.ll2004l : iBonanos et alJl2004 l. On the other hand, new statistical 



analyses indicate th at there appears to be an upper mass limit somewhere in t he 
120-200 M range dWeidner fc Kroupall2004t lOev fc Clarkell2005l : lFigedl200fift . 



In order to study the effect of unresolved binaries on the slope of the IMF at 
its upper-mass end and to obtain a better constraint on the stellar upper mass 
limit, we are currently engaged in an HST GO program (10602) which is a con- 
tinuation of a previous one (10205). We are obtaining multi-filter imaging of all 
known Galactic 02/03/03.5 stars with the High Resolution Channel (HRC) on 
the Advanced Camera for Surveys (ACS). The HRC has a pixel size of 0'.'027 and 
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Figure 6. HRC images of the core of Trumplcr 14. The left panel shows 
the F435W (B) image and the right panel shows the F850LP (z) image. The 
field size is 31.0" x 28.6". Top is 23.2 degrees West of North. Geometric 
distortion has been removed from the images but some cosmetic artifacts are 
still present in the F850LP case. 



its P SF is well sampled and very stable across the detector ( Anderson Sz Kingl 
2004). Furthermore, its geometric distortion is very well characterized, allowing 
for a relative astrometric precision of 0.005 pixels for very bright stars. 

We present here our first results on the core of Trumpler 14 (Fig. EJ, a 
young cluster in the Carina Nebula Association that contains at least three 
very-early O-type stars, including HP 93129 A, of spectral type 02 If *, the 
closest known 02 star ( Wal born etaD 120021 : [Mafz Apellaniz et al .1120041 ) . and 
HP 93129 B, which is an 03.5 V((f+)). Those two stars are sep arated by 2'.' 7 
and a ppeared to be single in ground-based s peckle interferometry ((Mason et al.l 
1998). Recently however. INelan et al.l (l2004h were able to split HP 93129 A into 
two components using the Fine Guidance Sensor (FGS) on HST. They obtained 
a separation of 55 ± 3 mas at a position angle of 356 ± 4 degrees (measured 
from N towards E) and a magnitude difference of 0.90 ± 0.05 in the visible. In 
the same data HP 93129 B is unresolved. 

We applied a PSF-fitting IPL photometry code especially written for this 
purpose to the HRC data for HP 93129 A and B. The code was applied to two 
dithered exposures in each of the F220W and F435W filters, thus yielding four 
independent measurements for each star. If a single component is used for the 
fit, the residuals for HP 93129 B are very small but those of HP 93129 A are 
very large, as expected for a binary system (Figs.[?Knd|Hl). On the other hand, a 
two-component fit yields very small residuals for HP 93129 A, hence confirming 
the binary character detected with FGS. 

We present the results of the two-component PSF fitting to HP 93129 A in 
Table El The final row shows the proposed values for the separation, position 
angle, and Am derived from the four independent HR C exposures. Our po sition 
angle is 18 ± 4 degrees to the E of that measured bv INelan et alJ ( 20041 ). sug- 
gesting that we are detecting the relative motion of HP 93129 Ab with respect 
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Figure 7. Postage stamps extracted from one of the two F435W HRC ex- 
posures for the two main components of HD 93129 (A, left and B, center) and 
PSF used for fitting (right). A logarithmic scale between 0.1% and 100% of 
the peak value is used in all cases. The field size is 0'.'53 x 0'.'53. Star symbols 
are used to identify the positions obtained by PSF fitting. 



Figure 8. Fit residuals for HD 93129 A (left and center) and HD 93129 
B (right) for one of the F435W exposures. The central panel assumes two 
components (e.g. it treats A as a binary system) while the other two assume 
a single component. A linear scale between —3% and 3% of the peak value 
for the data is used in all cases. The residuals have been smoothed with the 
PSF. The field size is 0'.'53 x 0'.'53. 



Table 3. PSF analysis for the two components of HD 93129 A. 



Separation (mas) Position angle (°) 



Am 



F220W exp. 1 
F220W exp. 2 
F435W exp. 1 
F435W exp. 2 
Final value 



50.4 ± 2.5 

51.4 ±2.5 

52.0 ±2.5 

52.1 ±2.5 

51.5 ± 1.2 



16.8 ± 2.8 1.1589 ±0.0087 

13.3 ±2.7 1.1751 ±0.0086 

14.4 ±2.7 1.1284 ±0.0055 
13.2 ±2.7 1.1244 ± 0.0052 
14.4 ± 1.4 1.1670 ± 0.0061 (F220W) 

1.1263 ± 0.0038 (F435W) 
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to Aa. Prom the 2.4 year difference between the epochs of the FGS and HRC 
observations we derive a very preliminary orbital period of ~50 years for the 
system but, obviously, observations at other epochs will be needed to confirm 
and measure an orbit. We point out that this is the first 02/03/03.5 star ever 
measured to be an astrometric binary. 

The magnitude difference between HD 93129Aa and Ab measured from the 
HRC data is similar to but slightly larger than the one measured with FGS 
(which uses light with longer wavelengths). The two components have very 
similar F220W-F435W colors, with Ab being redder only by 0.0407 ± 0.0072 
magnitudes. Given their proximity, it appears unlikely that the relative color 
is caused by differences in the amount or type of extinction. On the other 
hand, such a difference in color is equivalent to that between 50 000 K and 
44 000 K for TLUSTY models (lLanz fc Hubenvl I2003T ) with logg = 4.75 and 
solar metallicity, indicating that both components are likely to be early-0 stars. 
We are curre ntly working on a mo re detailed analysis of the photometry using 
CHORIZOS (|Mafz An ellanizl 12004 ) . 

It is important to note that HD 93129 A has not been identified as a spec- 
troscopic binary. With a separation between its two components of about 150 
AU, one would expect relative velocities of the order of 30 km s _1 if the inclina- 
tion is large. Therefore, its non-identification is not surprising, since one would 
require a separation an order of magn itude smaller to allow fo r a clear detection 
of radial velocity variations 1 (see e.g. [Bonanos & Stanek 2005). This also means 
that we cannot even discard the possibility that either HD 93129 Aa or Ab are 
binaries themselves. 

What does this mean for the biases in the top end of the IMF induced by 
multiplicity? Trumpler 14 is at an approximate distance of 2.7 kpc. If it were 
located at the same distance as the Galactic Center or NGC 3603, HD 93129 A 
will likely appear unresolved with HRC or FGS. At the distance of the Magellanic 
Clouds, Aa and Ab would be unresolved and A and B could be resolved but only 
with HST or adaptive optics. Moving to M31 or M33, HD 93128 (another 03 
star in the cluster outside the field in Fig. EJ) and HD 93129 would have an 
angular separation similar to that of Aa and Ab at its actual distance, with the 
rest of the stars (likely of late-0 and B type) in Fig. in between. This yields 
a total of (at least) four early-type stars blended together in a cluster that is 
quite massive, but far less than R136. The reader can easily deduce from these 
simple calculations how much he/she can trust IMF derivations of extragalactic 
young clusters that do not take into account multiplicity corrections. 



5. Conclusions 

We have discussed three sources of biases in the determination of the initial 
mass function (IMF): [1] the use of constant-size bins, [2] the uncertainty in the 
determinations of masses, and [3] the existence of unresolved multiple systems. 
Those three effects tend to produce IMFs that are flatter than the real one, all 



1 Oi course, close binaries are easier to detect not only because the radial velocity variations are 
larger but because they occur on shorter time scales. 
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of them large enough to potentially introduce significant systematic errors in 
the derived power-law slope. In the first case we present a technique that can 
get rid of the bias almost completely, even after the data have been obtained. 
In the second case we present a method that can also be used a posteriori to 
estimate the biases and we give some advice as to how to reduce them by using 
multifilter data. The third demonstrated by the example of HD 93129 A, 

is harder to correct, given that the current capabilities do not allow us to detect 
all multiple systems in the stellar clusters and associations of interest, even those 
in our own Galaxy. Nevertheless, that should not preclude us from making an 
effort to close the current gap between visual and spectroscopic binaries or from 
trying to estimate the contribution of of unresolved multiple systems to the IMF 
slope. 

Acknowledgments. Support for this work was provided by NASA through 
grants GO-09419.01, AR-09553.02, and GO-10205.01 from the Space Telescope 
Science Institute, Inc., under NASA contract NAS5-26555. 

References 

Anderson, J., & King, I. R. 2004, ACS Instrument Science Report 2004-15 

Bonanos, A. Z., & Stanek, K. Z. 2005, in ASP Conference Series, Vol. 332, 257 

Bonanos, A. Z., et al. 2004, ApJL, 611, L33 

Chabrier, G. 2003, PASP, 115, 763 

Dolphin, A. E. 2000, PASP, 112, 1383 

Fall, S. M., & Zhang, Q. 2001, ApJ, 561, 751 

Figer, D. F. 2005, Nature, 434, 192 

Kouwenhoven, M. B. N., Brown, A. G. A., Zinnccker, H., Kaper, L., & Portegies Zwart, 

S. F. 2005, A&A, 430, 137 
Kroupa, P. 2001, MNRAS, 322, 231 
Lanz, T., & Hubeny, I. 2003, ApJS, 146, 417 
Lucy, L. B. 2000, MNRAS, 318, 92 
Maeder, A., & Meynet, G. 2001, A&A, 373, 555 

Mafz Apellaniz, J., Walborn, N. R., Galuc, H. A., & Wei, L. H. 2004, ApJS, 151, 103 
Mai'z Apellaniz, J. 2004, PASP, 116, 859 

Mai'z Apellaniz, J., & Ubeda, L. 2005, ApJL (20 August 2005 issue) 

Mason, B. D., Gies, D. R., Hartkopf, W. I., Bagnuolo, W. G., Brummelaar, T. T., & 

McAlister, H. A. 1998, AJ, 115, 821 
Massey, P., & Hunter, D. A. 1998, ApJ, 493, 180 

Nelan, E. P., Walborn, N. R., Wallace, D. J., Moffat, A. F. J., Makidon, R. B., Gies, 

D. R., & Panagia, N. 2004, AJ, 128, 323 
Oey, M. S., & Clarke, C. J. 2005, ApJL, 620, L43 
Rauw, G., et al. 2004, A&A, 420, L9 
Stetson, P. B. 1987, PASP, 99, 191 

Ubeda, L., Mafz Apellaniz, J., & MacKenty, J. W. 2005, in preparation 
Walborn, N. R., et al. 2002, AJ, 123, 2754' 
Weidncr, C, & Kroupa, P. 2004, MNRAS, 348, 187 

Wheaton, W. A., Dunklee, A. L., Jacobsen, A. S., Ling, J. C, Mahoney, W. A., & 
Radocinski, R. G. 1995, ApJ, 438, 322 



